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Abstract 

The method of flow equations is applied to QED on the light front. Requiring that 
the partical number conserving terms in the Hamiltonian are considered to be diagonal 
and the other terms off-diagonal an effective Hamiltonian is obtained which reduces the 
positronium problem to a two-particle problem, since the particle number violating contri- 
butions are eliminated. No infrared divergences appear. The ultraviolet renormalization 
can be performed simultaneously. 
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1 Introduction 



We use the method of flow equations to construct an effective Hamiltonian starting from the 
light-front formulation which can be used to solve the bound state problem. 

Recently Glazek and Wilson QTJ and one of the authors have suggested methods to 
diagonalize Hamiltonians continuously which have been called similarity renormalization and 
flow equations by the authors, resp. 

It is common to both methods that they eliminate by means of a unitary transformation 
initially the off-diagonal matrix elements between states with large energy differences and con- 
tinue with states closer and closer in energy, so that off-diagonal matrix elements between states 
of energy difference larger than A are eliminated or strongly suppressed. The final aim is to 
eliminate them completely (A — > 0) and to obtain a diagonalized Hamiltonian. 

Application of this method to an n-orbital model has shown, that a literal use of this con- 
cept can lead to convergency problems [0. Instead one may leave the idea of diagonalizing 
immediately in favor of block-diagonalizing. If matrix-elements between states of equal parti- 
cle number are considered diagonal, then the procedure brings the Hamiltonian into a block- 
diagonal form, a procedure which worked much better in the case of the n-orbital model, where 
block-diagonalization with respect to the quasiparticle number (number of electrons above the 
Fermi edge plus number of holes below the Fermi edge) was performed. 

It becomes apparent from the calculations by Jones, Perry and Glazek |J on the basis of 
the similarity transformation, that this scheme works well down to energy differences of the 
order of Rydberg, but if one goes below, then contributions in higher orders in the coupling 
become important. 

Indeed it seems to be rather difficult to obtain bound states from plane waves by continuous 
unitary transformations. In eliminating only the terms which do not conserve the number of 
particles one postpones the diagonalization, but reduces the problem to one in the space of 
fixed particle number M. Thus for the positronium problem it is sufficient to determine the 
one- and two-particle contribution of the Hamiltonian for electrons and positrons. 

Basically the procedure is very similar to that of the elimination of the electron-phonon 
interaction |J which yields an effective attractive interaction between electrons responsible for 
superconductivity. Both the method of flow equations and the similarity remormalization M 
yield results different from Frohlich's original ones || but in very good agreement with more 
sophisticated methods. In QED it is the interaction of the electrons with the photons instead 
of the phonons which has to be eliminated. 

A basic advantage of the methods of similarity renormalization and flow equations in com- 
parison to conventional perturbation theory is, that one obtains normally less singular effective 
interactions. 



This procedure is similar to the Tamm-Duncoff Fock space truncation ||, [10], [II], |12| in the 
sense that also in this truncation particle number changing interactions are eliminated. 
We define an effective, renormalized Hamiltonian 

H eff = lim ff fl (A,A) (1) 

where Hb(X,A) is obtained by means of the unitary transformation U(X,A) from the bare 
Hamiltonian (regularized at the cutoff A) 

H B (X, A) = U(X, A)# B (A)C/t(A, A). (2) 

The unitary transformation U(X,A) is determined by the flow equations below. 
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Flow equations eliminate the far-off diagonal matrix elements in the 'energy space' [] (the 
matrix elements between the states with large energy jumps \Ei — Ej\ > A, where A is the 
'running' UV cutoff) but only for those blocks that change the number of quasiparticles. The 
value of A = A — > oo corresponds to the initial bare Hamiltonian, a finite A determines the 
effective Hamiltonian at an intermediate stage, for A = the elimination of the 'non-diagonal' 
part in the Fock representation is complete fig. (1). 

The method of flow equations in the perturbative theory frame is discussed in section 2. In 
section 3 we apply this scheme to QED on the light front. The effective Hamiltonian serves for 
the calculation of bound states, which due to the triviality of the vacuum state is now much 
simpler. 



2 Flow equations in the perturbative frame 

The flow equation reads 

d -^ = Hl\H{l)\. (3) 

This is the differential form of a continuous unitary transformation acting on the Hamiltonian 
in general; rj(l) is antihermitean generator of the transformation and I is the flow parameter. 
The aim of the flow equations is to bring the initial Hamiltonian matrix H(l = 0) to a diagonal 
(or block-diagonal) form. Finite values of the flow parameter I correspond to intermediate 
stages of diagonalization with the band-diagonal structure of Hamiltonian (similarity renor- 
malization) |[L3|| . The result of the procedure at I —> oo is a block-diagonal Hamiltonian. Also 
the transformation is designed in a way to avoid small energy denominators usually present in 
perturbation theory. 

What is the choice of the generator r](l) that performs such a transformation? Break the 
Hamiltonian into a 'diagonal' and a 'rest' part H = + H r . Then the prescription in is 

rj(l) = [H d (l),H r (l)]. (4) 

At this step the unitary transformation is defined. The only freedom left is the principle of 
separation into the 'diagonal' and the 'rest' part. It depends on the problem one wants to treat. 

Our goal is to transform the Hamiltonian into blocks with the same number of (quasi- 
particles. This means, that we define the 'diagonal' part H d as the part of the interaction 
which conserves the number of particles (electrons, positrons, photons), and the 'rest' H r as 
the particle number changing part. In the case of QED(QCD), where the electron-photon 
(quark-gluon) coupling is present, the number of photons (gluons) is conserved in each block 
of the final effective Hamiltonian. 

As a result of the unitary transformation new interactions are induced (see below). They 
are absent at I = and are generated as I increases. They also give rise to new terms in the 
generator of transformation r](l). This in its turn generates new interactions again. 

To be able to perform the calculations analytically we proceed in a perturbative frame and 
truncate the series assuming the coupling constant is small. 

As illustration of the method we consider QED on the light front. For any finite value of / 
one has 

H(l) = Hf + + H® + H® + ... (5) 



1 Under the 'energy space' we understand the basis of the free Hamiltonian corresponding to the sum 

of kinetic energies of single noninteracting particles. 
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where the superscript denotes the order in the bare coupling constant, ~ e n ; the indices 
'd' and 'r' indicate the diagonal and the rest parts correspondingly. The part is the free 
Hamiltonian, corresponding to the single particle energies with the structure in secondary quan- 
tization a^a, tfb, d^d, where a, b, d are the annihilation operators of the photons, electrons and 
positrons correspondingly; denotes the electron-photon coupling (of the type a^tfb); Hf 
is the second order diagonal part of the Hamiltonian, having the structure b^d^bd, tftfbb, d^d^dd 
(in the light front they correspond to the canonical instantaneous (seagull) and to newly gener- 



ated interactions in the diagonal sector in second order |14|, [115|1 ). Note, that the diagonal part 
in the flow equations is not only the free Hamiltonian but the full particle number conserving 
part of the effective Hamiltonian. The choice of only as the diagonal part gives rise to 
the band-diagonal structure of the effective Hamiltonian in each 'particle number' sector in the 



similarity renormalization scheme [I3|. However, this makes a difference for the diagonal part 
only if one goes beyond third order in e. 
The generator of the transformation is 

V (l) = [H d , H r ) = [Hi?\ «] + [Hi?\ H^} + ... = + r/ 2 ) + ... (6) 

Up to second order the flow equation reads 

= fo, H] = [[H$\ HU),H®) + [[Hf\ H(%m + [[tff, H(% H®) + ... (7) 



Also terms of higher orders in e are generated by the flow equations. 

AO' 



In the basis of the eigenfunctions of the free Hamiltonian 



Hf ] \i >= Ei\i > 

one obtains for the matrix-elements between the many-particle states 



Vij = ( E i - Ej)H^ + (Ei - Ej)H r ij + ... 

= -(Ei - EjfH^ + H^]ij - {Ei - E 3 ) 2 H { $ + ... (9) 



dl 

The energy differences are given by 

n.2 ni 



Ei — Ej = 22 E itk — ^ Ej,k (10) 
fc=i fc=i 

where E^ and Ej ^ are the energies of the created and annihilated particles, respectively. 

The energy Ei depends on the flow parameter I only in second order in the coupling. There- 
fore one has 

a "J = -(E- -E-?H {1) 
dl 13 

(I) = (I = o) e -(Ei-Ej)*i = H (X) (A = A ^ 00 ) e - 1£ ^ (11) 

Here we have used the physical meaning of the flow parameter /. Namely, in the similarity 
renormalization scheme it defines the width of the band A, corresponding to the UV-cutoff, 
where the matrix elements of the effective Hamiltonian are not zero (\Ei — Ej\ < A) [13[]. The 
connection between these two quantities is 

l = y 2 (12) 
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In the flow equations as used here the matrix elements of the interactions, which change the 
number of particles, are strongly suppressed, if the energy difference exceeds A, while for the 
particle number conserving part of the effective Hamiltonian the matrix elements with all energy 
differences are present. 

In the flow equations A characterizes the smooth UV-cutoff. This fact insures the analytical 
behavior of the effective Hamiltonian with A, that helps in numerical calculations. 

In second order one has to distinguish between the behavior of the 'diagonal' and the 'rest' 
term. For the 'rest' part one has 

d -^=[^\H^U-{E l -E ] fH^ 1 (13) 

where index V by [rj^\ H^'] r defines the particle number changing part of the commutator. 
Introduce 

^j(i)=e-^-^(§(0. (14) 
Then the solution reads 

H%{1) = H£](l = 0) + f Q dl'^- E ^ l \^\H^\ ij {l'). (15) 

For the 'diagonal' part one has 

J ^=[V (1 \H(%, (16) 
and the solution is 

= H^Al = 0) + f dl'[^\H^\ dij {l'). (17) 



L dij\' J ) dij' 

Note, that though in general the commutator [[H^\ Hjp], H^] is not zero, it is not present in 
the flow equation due to the definition of the diagonal part. The corresponding commutator 
\[H^\ H^]H d Q \ in the 'non- diagonal' sector insures the band-diagonal form for the 'rest' in- 
teraction and also gives rise to the different structure of the generated interaction (the integral 
term) in the 'rest' and 'diagonal' sectors. 

The commutator [r/^, H^] gives rise to new terms in second order in the bare coupling 
e. In the case of QED it induces new types of interactions and generates the renormalization 
group corrections to the electron (photon) masses. The coupling constant starts to run in third 
order in e. 

Note, that the method of flow equations (and also the similarity renormalization) enables 
one to build an effective low energy Hamiltonian together with all, 'canonical' and 'new' fl6| 



counterterms, found (from the coupling coherence condition [|T(J , |TJJ ) order by order in the 
coupling constant e. This defines an effective renormalized Hamiltonian which may be used for 
the numerical solution of the bound state problem. 

3 Renormalized effective electron-positron interaction 

In this section we give the effective Hamiltonian in the light front dynamics for the positronium 
system, generated by the unitary transformation in section 2, compare JTJ] . 
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The light front Schrodinger equation for the positronium model reads 



H LG \lp n >= M*|Vn > (18) 

where H^c — P^P^ is the invariant mass (squared) operator, refered for convenience to as the 
light front Hamiltonian of positronium and \if) n > being the corresponding eigenfunction; n 
labels all the quantum numbers of the state. 

The canonical Hamiltonian of the system Hlc contains infinitely many Fock sectors (i.e. 
one has for the positronium wave function \ip n >= c eS \(ee) n > +c eBj \(eej) n > +c eg77 |(ee77) ri > 
+...) and each Fock sector contains states with arbitrarily large energies. We now 

(1) introduce the bare cutoff (regularization) with the result H^ C (A) - the bare Hamiltonian; 

(2) perform the unitary transformation by means of flow equations eqs. (|1|) and (|2|) with 
the result H^fJ - the effective renormalized Hamiltonian (table 1 for finite value of A); 

(3) truncate the Fock space to the lowest Fock sector (\ee >) with the result H^fJ - the 
effective renormalized Hamiltonian acting in the electron-positron sector. 

Then the eigenvalue equation reads 

Hl f c f \(ee) n >=M 2 n \(ee) n > . (19) 

The effective light front Hamiltonian consists of the free (noninteracting) part and the effective 
electron-positron interaction 

H$ = flffi + V$ (20) 



The light front equation eq. ( |T9"| ) is then expressed by the integral equation (the coordinates 
are given in fig. (2) 

+ 53/ < x ' , k '±'i s 3, s 4 V lc\ x , k ±', s i, s 2 > ip n {x, k ± ;s 1 ,s 2 ) = (21) 

sT^2 JD 2 ( 27F ) 

The integration domain D is restricted by the covariant cutoff condition of Brodsky and Lepage 
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m 2 + k 2 , lL> , 



x(l — x) 



< A + Am 1 (22) 



which allows for states which have a kinetic energy below the cutoff A. 

For the effective electron-positron interaction one has in the exchange and annihilation 
channels 

vi f j = v exch +v ann = £ \imiyr+vr st +vn m 

channel 

where the terms in eq. (|2~3D correspond to generated, instantaneous and perturbative photon 
exchange interactions in the electron-positron sector. The interactions generated by the flow 



equations from the Hamiltonian in the light front frame |1J] are, fig. (2) 
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in the exchange channel 



yinst 



4e 2 



Ax 



r<5 



A2 s l s 3 



PT 



-e 2 iV 1 ^/ A (A 1 )/A(A 2 ) 
^3 



(24) 



in the annihilation channel 



yinst 



e 2 N 2 



V 



M 2 



+ 



Mo 2 



4e 2 5 - 5 - 



where /a (A) is the similarity function 

/a (A) = e"£ 



(25) 



(26) 



A 2 
P+- 



where we have used the rescaled value of the UV cutoff A 

As long as A is finite we have taken the still existing terms in into account by the 
perturbative result obtained by Lepage and Brodsky |18| . 

The functions Ni,N 2 (current-current terms) and the energy denominators Aj, i = 1,2,3, 
Mq, Mq 2 , M 2 are defined in the light front dynamics [15] as follows [Hj fig. (2) 



(x — x') 



i\2 



xx'(l — x)(l — x') 



+imv2(x' — 



x 



6siS3 6s2Si xx' T ' ' ^ + 6s ^ 6s2S * (1 _ x )(i _ X ') 

1 



No 



^•4 



+im\/2 



si 



xx'(l — x)(l — x') 

S3 



^S3S4^SiS2 „\^3 ' £ Si ^S3S4^SiS2 ,/ 1 A "^4 ' ^S4 



— Xj 



! x'(l — x') 



\/2 



>>*) 



(27) 



and 



2 (k ± -k> ± y | fcj( 82 ) | fc2(s 2 ) 



(x — x') (1 — x) (1 — x') 



X 



/u I (5) — I ~~r~ ISE'ijJv 



(1-x') 



r: 



4 (l-x) 



r 2 = 2 



fcj(£i) 
x 

s = — s 



(x — x') x x' 
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with the definitions 

(xk' ± — x'k^) 2 + m 2 (x — x') 2 



Ax 



XX 



A 2 — Ai| x _ > (i_ a! ) )a j_ > (i_ a! /) 



Ai 

1,7 



Ax 



A 9 



A 5 



x — x 



A 3 = _(A 1 + A 2 ) + |x-x'| -(M 2 + M, 



£; 2 i + m 2 



" " x(l -x) ' 
(^) 2 + M 2 
P+ 



2 + m 2 



,! ~~ x'(l - x') 
P= (P+,P ± ) 



M n = 2m + P n 



(28) 



Here x is the light front fraction of the electron momentum, P is the total momentum of 
positronium and the binding energy of the positronium. 
Explicitly the integrals in eqs. ( |2"4| ) and ( ^5P are given 



(29) 



We use instead of the light front parameterization fig. (2), the instant form. Namely, we express 
the variable (x, k±) in terms of the equal-time variable k = (k z , k±) as 



x 



1 + 



k 7 



k 2 + m? 



k — k 7 ~\~ k I 



(30) 
(31) 



and similarly for x' and k 2 as function of k z . 

Then the electron-positron interaction reads in the exchange and annihilation channels 



Vx 



-e 2 N hX 
+ 



'Ai + A 2 ' 
A? + A|, 



1-e aF*) + — e = ~\^~ 
A 3 



4e 2 



[x — x 



,\ 2 ^l s 3 ^S2S4 



+e 2 iV, 



'm 2 + m o 2 ' 



1 "o+v 



M$ + M' \ 

where one has in the instant frame 
M 2 = 4(A? 2 + m 2 ) 



M n+ M n 4 , 

1! 



(32) 



(33) 



and similarly for M' 2 as function of k 2 . For all the quantities, defined in eqs. ( p7|) and (|28|) , the 
substitution x(k z ), x'(k z ) is to be done. To the leading order of the nonrelativistic approximation 
\k\/m <C 1 one obtains (for the exchange channel) 



Ai ~ A 2 ~ A 3 = A = (k - k') 
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^" K - e2 ^V i_/l(A)) 

(34) 

X' — X 

This gives for the electron-positron interaction in the whole nonrelativistic range of A, A <C m 

V ~ 6 ( £_£)2 (x ™ 
the A- independent result. Making use of the following expressions 

N diag _ 4 ( k ± ~ k ±) 2 6 s 

iV l ^ (x — x') 2 

(£ - k'f = (k ± - k' ± f + (Jfe a - « (£ ± - ^) 2 + 4m 2 (x - x') 2 (36) 

one obtains in leading order of the nonrelativistic approximation the 3- dimensional Coulomb 
interaction (e 2 = Aira) 

l "" > M 16m ' ("(T^i) 4 " A "*' (37) 

hence the rotational invariance is restored in this order. This result eq. fl3"T| ) is valid for any 
nonrelativistic value of cutoff A and does not depend on the details of the similarity function 
/a(A). 

We perform the limit of A — > in the effective interaction, eqs. ( f24|) and (|25[) , that corre- 
sponds to the complete elimination of the electron-photon vertex. Then the perturbative term 
for the dynamical photon exchange V PT vanishes. Therefore the effective interaction eq. (p3"D , 
generated by the flow equations, is defined in the whole parameter region, (except maybe for 
the Coulomb singularity point q = k — k' = 0) as follows 



VlH — Vexch + V m 




Ai + A 2 \ , ( 4e 2 * 
A 2 + A 2 ,/ V (a-*') 2 1-8 , 

M 2 + ilC 
M 4 + Mo\ 



+ (4e 2 5 SlS3 8 S2 s 4 ) , (38) 



The effective interaction V^q obtained in eq. (pjl) stands for the kernel in the integral 
equation eq. (|2i~]) for the calculation of the bound state spectrum and the wave functions of 
positronium. The integral equation eq. ([H]) with the effective interaction given above is to be 
used for the numerical calculations of positronium spectrum. 

This interaction was also used for an analytical analysis of the positronium ground states. 
The standard singlet-triplet mass splitting for positronium ^a 2 Ryd was obtained and the de- 
generacy of the triplet ground state n — 1 was recovered [II], [EL 

Instead of the choice rj^ = (Ei — Ej)H^ one can choose = sign(Ei — Ej)H^. With this 
choice in the above equations the ratio nas ^° ^ e replaced by ? 2 A . If one approximates 
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M 2 = + M' q 2 ) in the expression for A3, then A 3 equals Al + A2 ; and one obtains the 

perturbation theoretic result V^T^^ + V mst |L8| for the effective interaction. These both 
interactions with different choices of i] do not coincide. Both have a leading Coulomb behavior, 
but they differ in the order e 2 q°, which determines the fine structure splitting. Quite generally 
the two particle interaction is of order q~ 2 and subleading terms in q. In the order of fine 
structure splitting a 4 also terms of order e 4 q~ 1 and e 6 q~ 2 will in general be important. Thus 
we expect that the difference in order e 2 q° will be compensated by differences in order e 4 and 
e 6 . They should yield in total the correct fine structure splitting. 

Using the similarity transformation Brisudova and Perry |19[ have obtained the correct spin- 
spin interaction for the positronium from the effective interaction in order e 2 . This had to be 
expected, since each spin enters the interaction with a factor of order q/m as compared to the 
leading Coulomb interaction. Thus the two-spin interaction enters only in order q° or higher. 
The only contribution to order a 4 comes from order e 2 . The same holds for the spin-triplet 
splitting (which is quadratic in the spin), and has been obtained correctly for the ground state 
by Jones, Perry and Glazek ||. By the same reasoning contributions to the spin-orbit coupling 
are of order q" 1 , so that contributions from order e 2 and e 4 have to be expected, but not from 
order e 6 . Spin independent contributions can be of order q~ 2 . To obtain all contributions of 
order a 4 one has to consider the interaction in order e 2 , e 4 and e 6 . 



4 Conclusions and outlook 

In this work we have outlined a strategy to derive an effective renormalized Hamiltonian by 
means of flow equations. Application of the flow equations with the condition, that particle 
number conserving terms are considered diagonal and those changing the particle number off- 
diagonal led as in other cases to a useful effective Hamiltonian. 

The main advantage of this procedure as compared with the similarity renormalization 
of Glazek and Wilson [|]] is, that finally states of different particle number are completely 
decoupled, since the particle number violating contributions are eliminated down to A = 0. 
Thus one is able to truncate the Fock space and the positronium problem reduces to a two 
particle problem which can be analyzed further analytically (since in leading order one obtains 
the nonrelativistic Coulomb problem) or numerically [?] for positronium bound states. 

The effective Hamiltonian, obtained by the similarity transformation, is band-diagonal in 
the energy space. The width of the band A introduces the artificial parameter in the procedure, 
which is defined from the physical reasoning (A is low enough to neglect the contribution of high 
Fock states, but is restricted from below to stay in perturbation theory region). Flow equations 
as used here with the particle number conserving part of Hamiltonian to be diagonal, have no 
additional parameter and converge well as A — > [F2J to the effective Hamiltonian, which is 
block-diagonal in particle number and can be used directly for the numerical calculations of 
the spectrum (work in preparation). 

The procedure of elimination of nondiagonal blocks, that change the number of quasiparti- 
cles, is performed not just in one step as in the method of Tamm-Dancoff truncation but rather 
continuously for the states with different energies in sequence. This is the main advantage of 
the proposed method as compared with Tamm-Dancoff truncation, the possibility to perform 
simultaneously the ultraviolet renormalization of the initial Hamiltonian. In general, in the 
definite order of perturbation theory all counterterms, associated with canonical operators of 
the theory and also with possible new operators induced by unitary transformation, can be ob- 
tained in the procedure [|13| . Since different sectors of the effective Hamiltonian are decoupled, 
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one does not encounter the usual difficulties of Tamm-Dancoff truncation and the methods 
related to it. Namely, the counterterms to be introduced are 'sector-' and 'state-' independent 
(work in preparation and |4]]). 

If one goes beyond the tree approximation then one obtains terms with ultraviolet diver- 
gences which have to be renormalized. This has not been considered in this paper. Simulta- 
neously also terms describing interactions between more than two particles are generated. In 
this approach we were not faced with infrared problems (except infrared collinear divergences 
along the light front, that are removed in the considered sector, if all possible diagrams in this 



order are taken into account [14 



By means of the flow equation method one can simultaneously renormalize the initial field 
theoretical Hamiltonian and construct the effective Hamiltonian, for which the Fock space trun- 
cation is valid. In order to solve the flow equations analytically we were forced to apply in this 
work the perturbative theory expansion. One is able to improve this approach systematically 
by going to higher orders in the coupling. 

We consider flow equations as a method which can also be used beyond perturbation theory 
in a self-consistent way. Examples in solid-state physics are the flow of the tunneling-frequency 
in the spin-boson model |7| and of the phonon energies in the electron-phonon coupling ||. 
Due to the flow the couplings decay even at resonance. 
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Figure 1: Flow equations perform the block-diagonalization of the bare Hamiltonian of the 
canonical theory ife(A) into a Hamiltonian consisting of blocks with equal number of particles. 
For a finite value of A the matrix elements of the 'particle number changing' sectors are squeezed 
into an energy band with roughly \Ei — Ej\ < A (left hand side picture) and are eliminated 
completely as A — * (right hand side picture). 

Figure 2: The effective electron-positron interaction in the exchange channel; the diagrams 
correspond to the generated and the instantaneous interactions. The perturbative photon 
exchange with two different time orderings is also depicted. 

Table 1: The effective light front QED Hamiltonian matrix up to second order in e in 
the Fock space representation. The matrix elements of the 'diagonal' (Fock state conserving) 
sectors are unrestricted in the energy differences; the 'rest' (Fock state changing) sectors are 
squeezed roughly in an energy band of width A. Black dots correspond to zero matrix elements 
in order 0(e 2 ). Instantaneous and disconnected diagrams are not included. 
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